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ABSTRACT 

Jupiter is expected to pulsate in a spectrum of acoustic modes and recent 
re-analysis of a spectrosco pic time series has i dentified a regular pattern in the 



spacing of the frequencies ( iGaulme et al.ll201lh . This exciting result can provide 



constraints on gross Jovian properties and warrants a more in-depth theoretical 
study of the seismic stru cture of Jupiter. Wit h current inst rumentation, such as 



the SYMPA instrument (jSchmider et al.ll2007l ) used for the lGaulme et al.l (l201l[ ) 



analysis, we assume that, at minimum, a set of global frequencies extending up 
to angular degree i = 25 could be observed. In order to identify which modes 
would best constrain models of Jupiter's interior and thus help motivate the next 
generation of observations, we explore the sensitivity of derived parameters to 
this mode set. Three different models of the Jovian interior are computed and 
the theoretical pulsation spectrum from these models for i < 25 is obtained. We 
compute sensitivity kernels and perform linear inversions to infer details of the 
expected discontinuities in the profiles in the Jovian interior. We find that the 
amplitude of the sound-speed jump of a few percent in the inner/outer envelope 
boundary seen in two of the applied models should be reasonably inferred with 
these particular modes. Near the core boundary where models predict large 
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density discontinuities, the location of such features can be accurately measured, 
while their amplitudes have more uncertainty. These results suggest that this 
mode set would be sufficient to infer the radial location and strength of expected 
discontinuities in Jupiter's interior, and place strong constraints on the core size 
and mass. We encourage new observations to detect these Jovian oscillations. 

Subject headings: Jupiter, interior 



Introduction 



Characterizing the interior structure of the Solar System's giant planets is critical for 
understanding their origin and evolution and the conditions in the protosolar cloud. Yet 
key characteristics of the interiors of the giant planets are still unresolved. We do not know 
the mass of Jupiter's and Saturn's core or whethe r the interiors of Uranus a nd Neptune are 
composed of compositionally distinct layers (e.g., 



Guillot and Gautierll2007f ). 



Measurements of their gravity fields can be used to discern the internal distribution of 
mass within the planets, however, these data are most useful in determining the density 
profile in only the outer portions of the planet. For Jupiter and Saturn, uncertainties in 
the equation of state for hydrogen at high pressures also substantially complicate efforts to 
disentangle the density profile within the inner half of the planet. For these reasons the mass 
of their pre sent cores is not well constrained. Possible masses range f r om 0-10 Earth masses 
for Ju piter ( Saumon and Guillot 2004 ; Fortnev and Nettelmann 201ol; Leconte and Chabrier 



2012 ). and 0-20 Earth masses for Saturn (lGuillotlll999l ; lGudkova and Zharkovlll999l ; ISaumon and Guillot 



20041 ). This lack of knowledge of core mass limits our understanding of the planetary forma- 
tion mechanism. The NASA Juno mission will in 2016 measure the gravitational harmonics 
of Jupiter to high accuracy, but these measurements will be most useful in constraining 
the structure of the outer regions of the planet, including for example, the nature of the 
deep atmospheric zonal flows. If Juno can measure the tidal response of Ju piter to Io, k 2 , 



the core mass might be measurable, but this is not assured (jHubbardl Il999l ) . Even if k 2 is 
measured, this parameter is highly degenerate with respect to central condensation. The 
same uncertainties that allow for different Jupiter models that all match J 2 exactly, will also 
apply to k 2 . In particular, internal layer bou ndaries introduce a degeneracy beyond that 
from uncertainties from the equation of state (jNettelmannl 1201 ll ) . Hence, an independent 
determination of layer boundaries in the deep interior is still required. 



With Jovian seismology, we have the unique opportunity to probe our interior modelling 
approach of assuming several homogeneous layers with sharp layer boundaries, in contrast 
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to an inhomogeneous, but smooth density gradient. If this approach should be invalidated 
for Jupiter, the hottest giant planet, we might conclude that it fails even more for the colder 
giant planets (Saturn, Uranus, Neptune) where, in addition, sedimentation is likely to occur. 

It has been recognized for well over thirty years that detection of aco ustic oscillations 



of a giant planet would provide a powerful probe of its interior structure (jVorontsov et al. 



19761 ). Since that time a number of papers have investigated the nature of Jupiter's os- 
cillation modes, and a number of attempts to detect Jovian oscillations have been made. 
With the recent det ection of the large global-mode spacing by ground-based observations 
( jGaulme et al.ll201l[ ). Jovian seismology has just begun to transition from its infancy to 
becoming a useful tool for constraining interior models. For comparison, helio- and stel- 
lar seismology have made tremendous strides in this same time period in both theory and 
observations, vastly improving our knowledge of the interiors of the Sun and some stars. 



Numerous auth o rs, for example 



Mossier et all jl99lh : [provost ; et all jl993h 



Vorontsov et al 



(119761) :IVorontsovl (ll98lh :lMarlevl (Il99lh: 



Gudkova and Zharkovj (Il999l ) , have attempted to 



study the frequency spectrum of Jupiter and/or Saturn both the oretically and observation- 
ally. On the theoretical s ide, the rapid rotation ( iVor ont sovl 1 1 9 8 lh and structure discontinu- 
ities ( iProvost et al.lll993l ) of Jupiter make for a challenging computation, but those modes 
which sense the core have frequencies of about 1 to 2 mHz, centered around 1.5 mHz. This 
frequency range depends fundamentally on the speed of sound in the metallic and molecular 
hydrogen envelope of Jupiter. Fortunately, the tools of stellar and solar seismology allow 
one to understand the best modes for constraining interior models. A better understanding 
of which particular set of frequencies will provide the greatest amount of information would 
allow more focused observational searches. 

Here we seek to compute an updated set of Jupiter oscillation frequencies for the most 
recent Jupiter models. We restrict the theoretical mode set to those with a spherical harmonic 
angular degree I up to 25 and frequencies below the acoustic-cutoff frequency. The choice for 
the maximum degree is roughly based on the spatial resolution of current instrumentation. 
Given this mode set, we then seek to establish what physical conditions in Jupiter's interior 
could be reasonably inferred using forward and inverse seismic techniques and given a few 
simplifying assumptions. 

The paper is organized as follows: The recently computed models are described in £j2j 
and their seismic properties are discussed in £j3j §H shows the computation of the forward 
problem, and a full description of the inverse problem, the solution that we choose to employ, 
and results of example inversions using the given modes is presented in §5J Finally, in §H1 
we summarize our main conclusions. 
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2. Jovian Models 

We have calculated the interior profiles of three different Jupiter models. The models 
are chosen to encompass the range of current knowledge. We assume the existence of a core, 
and explore core sizes up to the maximum currently allowed. Within the envelope, density 
discontinuities may occur, for instance, due to helium sedimentation. A discontinuity in the 
abundance of heavy elements may also occur at this location. The models are based on 
physical equations of state (EOS) for hydrogen and helium from different sources. In the 
following, we describe those models in detail. 



2.1. Model descriptions 

The SCV and LMR models 



M odels LMR and SCV are computed by the procedure described in iNettelmann et al. 



( 120081 . 120121 ) . They make use of the standard three-layer structure type, meaning two con- 
vective, homogeneous, and adiabatic envelopes above a core of heavy elements, where the 
abundances of helium and heavy elements in the hydrogen rich envelopes are chosen to 
reproduce observational constraints as follows. 

The He abundance in the outer envelope, which begins at the 1 bar level with tempera- 
ture Ti, is fixed to the Galileo entry probe value of Y% = 0.238(1 — Zx), where Z\ is the heavy 
element mass fraction there. In the inner envelope, the helium abundance, Y2, is adjusted 
to give an overall H-He mass fraction of 0.275 in agreement with the value of the protosolar 
cloud out of which the giant planets in the solar system formed. The abundances of elements 
that did not condense nor chemically bind to dust under protosolar cloud conditions should 
still be conserved in Jupiter, on average. The respective heavy element abundances in the 
envelopes, Z\ and Z%, are chosen to give an internal density profile that is consiste nt with 



the m easured gravitational moments J2 and J4 within the observational error bars (IGuillot 



2005|). 



While the position of the transition between the envelopes, parametrized by a transition 
pressure P1-2 influences the resulting solutions, especially the core mass M c , it is actually 
not fixed by an observable. In particular, we have P1-2 = 2 Mbar, M c = 1.28M®, J4/IO -4 = 
-5.84, Tx = 165 K, Y 2 = 0.2593, Z 1 = 0.0793 and Z 2 = 0.0817 for SCV, and Pi_ 2 = 8 Mbar, 
M c = 3.57 M ffi , J 4 /10" 4 = -5.89, T x = 170 K, Y 2 = 0.31 12, Z^ = 0.0381, and Z* = 0.1283 



for LMR. The LMR model is identical to model Jll-8a of INettelmann et all (12012 ) . 



Most importantly, models SCV and LMR differ in the equations of state (EOS) used 
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for H, He, and metals in the envelopes. Fo r model SCV we use d the SCvH EOS with H- 



SCvH-i for hydrogen, He-SCvH for helium ( jSaumon et al.lll995[ ). and this He EOS scaled 
in dens ity by a factor of four to repr esent heavy elements. For model LMR we used LM- 
REOS (INettefmann et al.ll2008l . |2012| ) with water represe nting metals. The pre s sure-d ensity 
relation in the cores was obtained from the rock EOS by lHubbard and Marleyl (119891 ) . 



To get smooth derivatives dP/dp with model LMR, the planetary P(p) profile in each 
layer was reduced to a small number (~ 10 — 40) of points Pi, transformed to an equidis- 
tant grid in logP on which then the derivatives (dlogp/dlogP)p i were taken by the 4-point 
Richardson method, and related to the radial coordinate values by spline or linear interpo- 
lation so that finally, c s = \J (P/ p) dlogp/dlogP, could be calculated. 



The MH08 model 



Model MH08 is the Jupiter model of iMilitzer et al.l (120081 ). with a 16M e core of ices and 
rocks, Y\ = 0.238 and Z 1 = 0.014 throughout the single envelope, and J 4 /10~ 4 = —6.14. The 
numerical derivatives along the published tabulated profile turned out sufficiently smooth. 



Summary of models 

To our best knowledge, these three models represent the widest range of currently 
offered Jupiter models with physical EOS. They span the observational uncertainty in the 
1 bar temperature (165—170 K); differ in core composition (pure rocks or ice-rock); include 
models with a small core (LMR) and a large core (MH08); have a strong, small, or absent 
layer boundary in the envelope; assume a helium enriched deep envelope (SCV, LMR) or 
not (MH08); assume that differential rotation has a negligible effect on J4 (SCV, LMR) or 
a significant effect (MH08), have atmospheric metallicities in agreement with the measured 
abundances (MH08) or with the prediction that oxygen is ~ 2 — 4 times solar (SCV, LMR); 
and, finally, include H-He mixing effects in the EOS (MH08) or not (SCV, LMR). 



2.2. Density and sound-speed profiles 

Figure [TJ shows radial profiles of the sound speed and density of the three models. In 
homogeneous layers, the sound speed c s (r) rises steadily with depth as the interior becomes 
warmer and increasingly degenerate. When dissociation occurs, the heat capacity at constant 
volume, Cy rises more than Cp, the heat capacity at constant pressure, and the rise of c s (r) 
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0.2 0.4 0.6 0.8 1 

Fig. 1. — (a) Density and (b) sound-speed profiles for the three models as a function of 
scaled Jupiter radius, Rj. Noted by arrows are small jumps in density for the LMR model 
at ~ 0.63i?j and for the SCV model at ~ O.SRj. Similar features are also present in the 
sound speed. 
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with depth is delayed. This effect can be seen at around 0.9 Rj in the MH08 and LMR 
models. The inner/outer envelope transitions of the LMR and SCV models are denoted by 
arrows in Fig. (U and will be the targets of our inversion analysis in the following sections. At 
the deepest layer boundaries, the steep rise in mass density maps onto a steep decrease in c s . 
The lower sound velocities in the deep interior of Jupiter as predicted by models SCV and 
LMR compared to MH08 are a consequence of higher envelope metallicities, hence higher 
mass densities there. 



3. Model seismic parameters 



Model profiles permit the computation of the expected oscillations in Jupiter, which 
are likely acoustic in nature and driven by convection, similar to the case of the Sun and 
solar-like stars. The oscillations are expected to couple to the troposph ere and produce 



small-amplitude albedo variations ( lMosserlll995t iGaulme and Mosserll2005l ) that in principle 



could be detected with appropriate observations. In what follows, we ignore rotational effects 

whose effects have been studied previously (e.g., 



Vorontsov 


1981; 


Mosser 


1990; 


Lee 


1993) 



frequency spectrum for non-radial pulsations (on the order of 27/xHz), our goal here is to 
provide an overall study on the usefulness of the pulsations on structure determination and 
these small deviations will not affect our main conclusions. 



3.1. Adiabatic oscillations 



We consider a linear theory of adiabatic, radial and non-radial oscillations. Compu- 
tations are carried out using small perturbations to the ord inary differential e quations of 
stellar structure, following the general formulation found in lUnno et al.l (119891 ) . For each 
model, we obtain the eigenfunctions and eigenfrequencies u® e for any given angular degree 
£ and radial order n. The superscript "0" distinguishes model frequencies from observed 
ones, as detailed later. In this study, we truncate the algorithm at I < 25. This angular 
degree limit is consistent with what c ould be obtained, in p rinciple, with the recently built 
and employed SYMPA spectrograph (ISchmider et al.l 120071 ) used for several ground-based 
observing campaigns. Of course, higher-degree modes should certainly be possible to detect 
from space on appropriate instruments with sufficient spatial resolution of the Jovian disk. 
For comparison, the Sun reveals degrees up to several thousand in high-resolution images. 



Fi gure [5] shows the theoretically computed cyclic frequency 



cu® e /2Tr as a function 
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Fig. 2. — Computed frequencies v® e for each of the three models in i — v space for I < 25 
and v < 3.5 mHz. Several radial orders (n) are labeled for reference. The dashed line along 
the lowest ridge denotes the dispersion relation for the surface-gravity (fundamental n — 0) 
mode. Black squares highlight the four modes for which kernels are plotted in Fig. 
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of the angular degree I for each of the three models. The lowest frequency "ridge" is the 
surface-gravity (fundamental) mode, which is not an acoustic mode and has a dispersion 
relation oj = yfgk, where u is the angular frequency, g is surface gravitational acceleration, 
and k is wavenumber: k 2 Rj = t{t + 1). The fundamental mode is analogous to deep ocean 



waves (jLighthilllll978l ). The rest of the spectrum shows the acoustic pressure modes up to 
I = 25. The maximum frequency we consider is 3.5 mHz, which is chosen as a value slightl y 
higher than the expected atmospheric acoustic cutoff frequency for Jupiter (]Mosserlll995l ). 
Physically, modes with higher angular degree are more sensitive to the outer envelope of 
the planet, while those with low degree probe the core. Increasing frequencies for a given 
angular degree denote successive radial orders n. The highest frequency modes correspond 
to roughly n w 22. The I = modes are the radial solutions, whose eigenfunctions have no 
nodes on the surface. 

It is unknown whether all modes that we expect to exist in Jupiter based on interior 
models will in fact be excited, or excited strongly enough to reach detectable levels. Connec- 
tive mode-excitation mechanisms are complex and a full understanding of all the processes 
responsible, including damping, is still lacking. 

We finally note that the solution finder of the differential equations is sensitive to the 
input frequency guess, and a solution is not always found for a given i since the input 
grid needs to be sufficiently fine. Thus, for different models there are occasionally "missed" 
modes, as one notices in Fig. As the equations are adiabatic, we also cannot provide the 
expected amplitudes of any p ulsation. We point out, however, that the recent observational 
work by iGaulme et al.l (120111 ) found that the highest amplitude mode had a frequency of 
about v = 1.2 mHz. This seems plausible, as recent work for the Sun and solar-like stars 
has shown a robust relationship b etween the frequency at maximum power (so-called ^ max ) 
and the acoustic cutoff frequency ( iBelkacem et al.l 1201 II ). In the case of Jupiter, z/ max should 
there fore be about one half of the expected acoustic cutoff frequency of v ~ 3 mHz ( jMosser 
1995h . 



3.2. Large frequency separation 

Various diagnostics computed from the seismic frequency spectrum are very useful to 
understand global properties that could be obtained from observations. To leading order for 
low-degree modes of high radial order, acoustic oscillations in convectively unstable bodies 
can be shown to follow a regular pattern given by 

Vnt— (n + - + i + a j Av, (1) 
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where a is a function capturing effects of near-surface regions (see lAerts et al.l l2010l and 
references therein). The quantity Av is the large frequency separation, and is the inverse of 
twice the sound travel time from the planet center to the surface: 



Av 



dr 



(2) 



The large frequency spacing quantifies the expected asymptotic spacing (for large radial 
orders n, i.e., high frequencies) between pairs of modes (n, £) and (n + 1, £), and scales with 
the square root of the mean density of the planet. The regular pattern in the asymptotic 
limit is somewhat visible by eye in the frequency spectrum in Fig. |2j 

The large frequency separation for the computed modes, along with the expected theo- 
retical value from Eq. [2] are shown in Fig. El For model LMR, we find Av = 156.1 /iHz, for 
MH08, Av = 161.7/aHz, and for SCV, Av = 155.5 /iHz. Asymptotic scaling relations for the 
Sun and solar-like s tars point to a large frequency separation Av = Av Q (M/M )°- 5 (i?/i? Q )~ L5 
(IBrown et al.lll99ll ). Given that Jupiter is about 1000 times less massive and about 10 times 
smaller in radius, one might expect the large frequency separati on for Jupiter to be r oughly 
similar to that of the Sun: Av & ~ 135 /aHz. The observations by lGaulme et al.l (120111 ) found 
Av pa 155.3 ± 2.1 /iHz for Jupiter, indeed very close to the three models studied here. 

Due to the equal spacing between modes of consecutive radial orders of the same de- 
gree predicted by asymptotic theory, vertical alignment of modes is expected if plotted in 
an echelle diagram, as is commonly done using frequency observations of stars. Figure H] 
shows such diagrams, where the frequency spectrum is vertically stacked from sections of 
a given length in frequency determined primarily by Av. As can be seen, Jovian models 
do not generally dem onstrate this vertical behavior for high frequencies (as first noted by 
Vorontsov et al.lll989l ). likely due to very strong core contributions to the eigenfunctions from 
the low-degree modes and sharp density jumps. The overall features of the echelle diagrams, 
particularly the wiggles and the flattening at low frequencies, are similar to wha t has been 
seen in previous Jovian models (jProvost et al.lll993l : iGudkova and Zharkovlll999l ). 



3.3. Small frequency separation 



From Eq. [T] it is apparent that frequencies of pairs of modes of (n, £) and (n — 1, £ + 2) 
should be similar in value. Expanding to second order in the a symptotic low-degree regime 
provides a correction to the mode frequency v n £ (jTassoullll980l ). This results in a breaking 
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f»+M ~ v "-< [ mHz ] v n +\,t - v n ,t [mHz] v„+ij - v„j [mHz] 

Fig. 3. — Large frequency separation Az/ of low-degree modes for each of the three models. 
The vertical dashed line is the theoretical value of Az/ for a given model from Eq. |2j 




i/^/Af - n - int(£/2) f llf /Af - n - int(£/2) iv/Ai/ - n - int(£/2) 

Fig. 4. — Echelle diagrams from the computed linear adiabatic modes for each of the three 
models. The symbols represent the different (low) angular degree of the modes. These plots 
can be thought of as stacked sections of roughly length Az/ of the frequency spectrum. 
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of the degeneracy of these pairs, and gives what is known as the small frequency separation 

n Az/ f Rl dc s dr 
4vr 2 z/„ £ J dr r 

This combination isolates contributions from the sound speed gradient in the integrand, 
which is most pronounced at the core-mantle boundary for the case of Jupiter, and therefore 
provides a strong diagnostic for that region. The pairs of modes we consider are the lowest 
degrees £ = (0, 2) and £ = (1, 3) for consecutive orders. The computed separations from the 
model frequency spectra and the expected separations from Eq. [3] are compared in Fig. 
The model separation is seen to span a range of about 100 /iHz and it is again evident that 
Jupiter does not simply follow the expected behavior predicted by asymptotic theory. The 
small frequency separation is also on the order of rotational effects on non-radial modes. 



3.4. Probe depths of the acoustic modes in Jupiter 

Modes excited in the convective envelope probe into the planetary interior to depths 
based on frequency and wavenumber, and refract back towards the surface due to the in- 
creasing density. These trajectories can be thought of as simple ray paths. To guide the 
inversion input, we show in Fig. [6] the lower turning point r t for modes of varying angular 
degree and frequency for all three models. The lower turning point can be understood as the 
region above which a particular mode is trapped beneath the surface. It is found through 
the relation 

C 1M = u2 u) 

r\ £(£+lY 1 } 

where the angular frequency u = 1i\v. Equation H] shows that, for a given £, higher-frequency 
modes probe deeper than low frequency ones, as Fig. [6] confirms. High-degree modes are 
trapped in the outer layers. In the analysis that follows, we exploit the basic property of 
seismology: that different modes sense different depths in the interior of the planet, to do 
targeted studies of distinct interior regions. 



4. Forward modeling - Sensitivity kernels for Jupiter 

Jovian oscillation frequencies are related to internal properties such as sound speed and 
density structure in a rather complicated and non-linear fashion. As is done in the cases of 
terrestrial and solar seismology to approximate such a relationship, we perform a linearization 
about an initial reference model, or three models, as is the situation here (for full details in 
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the solar context, see IChristensen-Dalsgaardll2002l ). For example, let po and Co denote the 
background density and sound-speed profiles of one of the models. Such a model permits 
calculation of a spectrum of adiabatic oscillations frequencies uj® t , as already depicted in 
Fig. [2J Now let us assume for the moment that several observations of Jovian frequencies 
are available, wj. We are interested in seeking corrections to the density, 5p = p — po, and 
squared sound speed, 5c 2 
modeled frequencies, Sco n £ 
linear relation such as 



0: 



UJ 



obs 
ni 



that match the differences between the observed and 
Assuming small differences to those quantities, a 



UJ 



5u, 



Writ 



R 3 r 



2 5c 2 



+ Kf(r) 6 Ar) 



dr + e 



(5) 



can be found. The quantities K are known as sensitivity functions, or kernels, and are com- 
puted using the model profiles and mode eigenfunctions. They are derived by straightforward 
first-order perturbation theory applie d to the basic equations of hyd rodynamics. Details and 
explicit expressions can be found in iGough and Thompson! (1199 lh . The sensitivity kernels 
relate the interior conditions to the frequency shifts Su. The first kernel above is for squared 
sound-speed differences (5c 2 ) without density differ ences, and the second is for d ensity dif- 



ferences (Sp) without sound-speed differences (e.g.. IChristensen-Dalsgaardl 120031 ). Relative 



uncertainties in each observed frequency due to realization noise are captured by e n £. 

Several example kernels for each of the three models are given in Fig. [7J Kernels for 
squared sound-speed are always positive, while those for density oscillate in sign. While 
overall these functions for each of the models are all quite similar, the LMR kernels for low 
angular degree have strong core-boundary sensitivity. In general, kernels for high-degree 
modes are less sensitive to inner regions, as expected from studying their lower turning 
points (Fig. [6]). We stress that the kernels are more useful than simply studying mode 
eigenfunctions, since they are computed for the specific observations in consideration: for 
this case these are frequency differences. 



5. Inversions of Jupiter frequencies 

Assuming there are several available frequency measurements, the associated uncertain- 
ties, and a set of computed kernels, Eq. [5] has just two unknown quantities, Sc 2 /c 2 and 
Sp/p. This is an inverse problem, that could in principle be solved by classical optimization 
methods, e.g., a least-squares algorithm. However, we emphasize that even in the absence 
of observations, as is the case here, detailed analysis can be done upon "inverting" Eq. [5] to 
determine important predictive properties, such as the internal resolving power and expected 
errors of the final inferred quantities. 
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We approach this problem by u tilizing a Subtractive Optim ally Localized Averages 
(SOLA) inversion procedure (see e.g. iPijpers and Thompson! fl992h . first developed for the 
case of the Sun. As shown below, it is a very useful tool even without available observations. 
For completeness, we fully describe a one-dimensional version of this inversion procedure 
that is constructed to give the most information at a given depth in Jupiter's interior. We 
follow with inversion examples using the Jovian models. 



5.1. Formulation of the SOLA technique 

Qualitatively, the main goal of the SOLA inversion is to combine the kernels in such 
a way that a particular region of the Jovian interior is isolated. As is evident in Fig. 
each kernel has a different depth sensitivity based on its frequency and angular degree. One 
can imagine that by adding and subtracting various weighted kernels, a desired averaging 
kernel can be constructed that peaks at one particular depth, has tails that tend to zero, 
and is sufficiently narrow to provide a good spatial resolution. Ideally, a delta function 
would be preferred; however a finite set of oscillation modes and the presence of noise makes 
that impossible. The coefficients that are used for this optimized linear combination of 
individual kernels are then use to combine the relative frequency differences in the same 
manner, ultimately providing an estimate of the unknown interior quantity. The SOLA 
inversion finds such a set of optimal coefficients. 

Assume M number of perturbations p (e.g., in this work M = 2, for p 1 = Sc 2 /c 2 
and p 2 = Sp/p) and N observations (i.e., N frequency measurements of identified modes). 
Equation can be rewritten as 

/ VA7(r)p 7 (r)dr + e h for i= 1...N, (6) 
Jo ~-i 



7=1 

where the subscript % corresponds to all n£ pairs of modes. We are interested in inferring a 
particular perturbation p a at some radial depth r , and according to SOLA 

P>o) = f>?(r ) — , (7) 

i=i 

where p a is an estimate of the true p a , and a a is some set of coefficients yet to be determined. 
In other words, all available measurements are averaged in some weighted fashion to infer 
the desired perturbation (a) at a single depth (r ). Plugging Eq. [6] into Eq. [7J gives 

N Rl M N 

p a (r ) = J>?(r ) / ^#7(rMr)dr + 5>?(r K 

i=l ^° 7=1 i=l 
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= / A a > a (r,r )p a (r)dr + / V A a ^(r, r ) p\r) dr + V <(r K (8) 
where we have defined the averaging kernel 

N 

A^(r;r ) = Y,at(r )K;'(r). (9) 

i=l 

Equation [8] is written in such a way as to reveal an important point: the estimate p a of 
the real p a from the inversion is given by three terms. First, it is the convolution of the 
true unknown quantity with the averaging kernel. If the averaging kernel is a well localized 
function centered around r , then this is a desirable result. The second contribution is from 
the contamination of the other perturbations 7 7^ a to our estimation, known as "cross talk," 
and a large value is not desirable. Ideally, one would like A a,J (r) = 0. The third contribution 
is from the propagation of noise through the inversion. 

To obtain the coefficient set a a , we minimize a cost function defined as 

Rj M N 

X °(a0= / r^lriroj-r^rj^f + ^a^roJA^^ro), (10) 
Jo 7=1 ij 

where T is a "target" function whose form is often taken to be 

2d 2 



T°"(r; r ) = exp I ^- ) 5« 7 . (11) 



It is thus a Gaussian and peaked at a target depth r = ro, where d controls the width of the 
function. It is only nonzero for the perturbation being sought, as denoted by the Kronecker 
delta function 5. The quantity A^- in Eq. [10] denotes the noise covariance matrix of the 
measurements. 

Minimization of Eq. [10] can therefore be interpreted as a competition between an av- 
eraging kernel that matches a prescribed target function and the contribution from noise. 
The quantity fi is a trade-off (free) parameter that allows one to tune between these two 
situations, depending on the type of solution desired. Ensuring a small "misfit" between the 
averaging kernel and the target can help localize around a particular depth, but will result 
in large error magnification, while relaxing that requirement gives more uncertainty and a 
coarser resolution, but a smoother solution. Examples will follow in the next section. 

Finally, to allay the detrimental cross-talk contribution seen in Eq. [HI we typically impose 
the constraint that 

Rj 

A a "(r;r )dr = 5 a7 . (12) 
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Therefore, at least the total spatial integral of the undesired averaging kernels vanishes even 
though it typically does not vanish at each radial point r. A second constraint is one that 
ensures the mass of Jupiter and the mass in the model are equal: 

4vr f R ' ^l p ( r y dr = . (13) 

Jo P( r ) 

The minimization of Eq. [10] with these two constraints leads to a system of iV + 2 linear 
equations that are solved using well-known matrix methods. 



5.2. Inversion examples 



Despite the very promising observational analysis in iGaulme et al.l (120111 ) in recovering 
the large frequency separation, no measurements with unique Jovian mode identifications 
are yet available. We proceed with example inversions using theoretical sets of modes from 
each model for the ranges i < 25 and < v < 3.5 mHz. The goal is to determine which 
observations would be most useful for inferring the sound speed and density profiles around 
several of the important interior regions of Jupiter, and by quantifying the spatial resolution, 
noise estimation, and overall internal sensitivity expected from such modes. The analysis and 
conclusions assume that the given models are reasonably close to the true Jovian interior. 

An estimation of measurement uncertainties is needed for the input to the inversion, 
and must therefore be approximated. We assume for simplicity that measurements are 
uncorrelated and have a relative uncertainty ti = 1%. A diagonal noise covariance matrix 
of identical elements a 2 will result, and the last term in Eq. [10] reduces to At of a?. The 
choice of this conservative value is based on the fact that at least such precision is required to 
identify Jovian modes so inversions can indeed be carried o ut. It corresponds ap proximately 
to Az//10, and is consistent with uncertainties provided in IGaulme et al.l ( 1201 ll ). 



5.2.1. Targeting the outer- envelope transitions 

We first consider the outer layers of the models where the envelope transition occurs. 
For models SCV and LMR, the radial value of the transitions are at ~ 0.8i?j and ~ 0.63-Rj, 
respectively. Model MH08 has only a single envelope and will be used subsequently for 
studies of the core. 

The jump in the SCV model sound-speed profile at ~ 0.8i?j is about 2.3% (see Fig. [Q. 
Figure [8] shows inversions for sound speed targeting this transition region. Panel (a) shows 



-17- 



the comparision of the two averaging kernels and the target function, centered at depth 
r = 0.8-Rj, with a full- width at half-maximum (FWHM) of 0.2i?j. The averaging kernel 
for the sound speed matches the target function very well, while the density perturbation 
averaging kernel is nearly zero throughout, thus eliminating any contributions of cross talk 
described earlier. Cancellation of sound-speed sensitivity is achieved in deep layers, while a 
small near-surface component remains. As a reminder, in all inversions the integral of the 
averaging kernel for sound speed is unity, while the integral of the density averaging kernel 
is zero, as constrained by Eq. [T2J 

The trade-off parameter in the inversion in Fig. [B^i (and all other inversions) is chosen 
so that a good match with low noise is obtain ed. This is done by in specting an "1-curve" 



for the point of optimal error and misfit values (IJackiewicz et al.ll2007l ). The combined noise 
after inversion is 0.01%. So even though we assume a 1% noise level in each measurement, 
the optimized combination found in the inversion is much lower, as one would expect. 

The magnitude of the inversion coefficients Oi(ro) is shown in Fig. [Bb as a function 
of lower turning point r t to help visualize which modes contribute most to this particular 
inversion. We see that those whose r t lies near the target depth seem to have the largest 
amplitudes. We also plot the magnitude of the coefficients in an I — v space representation. 
It is interesting to note that for some radial orders all modes along the ridge contribute with 
the same sign, while for other orders there are both positive and negative coefficients. 

The second inversion shown in Fig. $p and Eli is for a target function 4 times narrower, 
0.05i?j, with everything else fixed. To achieve a reasonable match to the target function, 
however, we find the noise increases to about 8 times the value found for the inversion with 
the wider target. This is more evident in the plotted magnitudes of the inversion coefficients, 
and demonstrates that it requires more extreme weighting of the kernels to find a match to 
the narrow target. Even such, the estimate of the sound speed at this target depth would 
have a standard deviation of only a ~ 0.04% (note that the noise values provided in the 
figures are for the squared sound-speed differences (c 2 — Cq)/c 2 : the true fractional error is 
approximately found by halving the number), making a detection of a 2.3% sound-speed 
jump as predicted by model SCV a realistic, although not a certain, goal. However, to 
reduce the noise to lower value, a trade off with the the spatial resolution would be need to 
be made, and more significant contributions from A p would factor into the inferences. 

The conclusion of this example is that detecting the small 2.3% jump in sound speed 
is reasonable to expect with this mode set based on noise levels, and the radial resolution 
appears to be sufficient to pinpoint such a narrow feature. The example also demonstrates 
the rather subjective balance between spatial resolution and noise that must be carefully 
chosen. 
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Figure [9] shows example inversions targeting the ~ 5.5% sound-speed jump in model 
LMR at ~ 0.63i?j. The top panels present the results for a wide target function of 0.2i?j. 
The sound-speed averaging kernel retains a sharp feature at this depth - there are not enough 
modes to "average" it away. The noise is acceptable at 0.004%, and the cross talk between 
sound speed and density again is minimal. Note the mostly positive sign of the coefficients 
for modes with lower turning point just beneath the target depth, and negative from those 
just above. It is thus critical to have a full mode set to get this type of localization. 

For the inversion in Fig. [Ufc and EH with a narrower target and therefore better overall 
spatial resolution, the noise is magnified by an order of magnitude from the previous case. 
For a sound-speed jump this large (~ 5.5%), however, the uncertainty (0.2%) would not 
necessarily mask the detection. Given all of our assumptions, we conclude that this would 
be a detectable signal. 

To emphasize how important it is to have as large a set of observed modes as possible, 
we repeat the inversion discussed in Fig. [9fc with a more limited mode set. We constrain 
the input to modes of £ < 10 and with frequencies between 1 and 3 mHz. The resulting 
inversion is shown in Fig. [TDJ The degradation in the potential inference is striking, as the 
lack of a sufficient number of high-degree modes greatly limits the inversion accuracy at this 
particular depth. 

5.2.2. Targeting the core-boundary transition 

We shift the focus to the core-envelope boundary and study the massive-core (16M e ) 
MH08 model. From Fig. [T] there is an overall ~ 80% drop in the density from the core region 
to the envelope. While this should produce a strong seismic signal, there are indeed fewer 
modes that probe to such depths. 

Figure [TT] shows inversions for density in model MH08 around r = 0.1i?j. We compare 
a series of inversions in which the FWHM of the target function decreases from top to bottom 
panels. We verify from the plots of the coefficients that the low-degree modes contribute 
most to the averaging at this depth. It is precisely because of the lack of a large number of 
low-degree modes that causes some difficulty in resolving the deep interior regions. There 
are also small but possibly significant contributions from cross talk with the sound speed, 
shown by the oscillating features in the deep interior. 

We next perform inversions after choosing a fixed target function with FWHM = 0.1i?j, 
and sweep through several values of the trade-off parameter \i to study its effect on the 
inversion. Figure IT21 shows the corresponding inversions in this case. The value of /i decreases 
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by two orders of magnitude in each panel from top to bottom - large /x means more emphasis 
on obtaining low noise, while smaller values typically gives better fits of the averaging kernel 
to the target, as Eq. [TD] describes. This trend is evident in the averaging kernels for density 
as the target/kernel match improves in each successive panel, while the noise increases as 
well. 

The maximum noise reaches about 1.2%, which is larger than the value of each individual 
measurement (1%). However, the spatial resolution is quite high and the cross talk is minimal 
for that particular inversion (even though it is clear from the amplitudes of the weights that 
most of the modes are weighted very strongly compared to the other examples). For a deep 
density jump of 80%, we would therefore expect to detect its location with this small set of 
modes, even though these uncertainties might not allow a fully accurate determination of its 
overall amplitude. 

We finally note that in these deep-layer inversions the high-degree modes are useful 
even though they do not probe deeply and only contribute rather small- amplitude weighting 
as Fig. [T2] shows. We have tested similar inversions using only £ < 10, and find that the 
fit worsens due to the fact that near-surface cancellation does not occur as efficiently than 
with these modes and non-local effects away from ro are significant. Overall, it is likely 
advantageous to use as many modes as possible in each inversion, regardless of target depth. 

6. Conclusions 

Recent seismic observations and analyses have indicated that the detection of pulsation 
modes of Jupiter is now possible, and the large frequency separation between low-degree 
modes of consecutive radial order has been measured. In this work we study three recent 
state-of-the-art models of the Jovian interior that span a broad range of parameter space. 
We assume a realistic acoustic oscillation mode set that, in principle, could be obtained from 
current ground or space-based instrumentation. This set contains angular degrees I < 25 and 
frequencies below the approximate acoustic-cutoff frequency of Jupiter of v < 3.5 mHz. We 
compute these modes using well-known methods, and find that the theoretical frequencies 
are consistent with past studies. 

For the first time, we show theoretical sensitivity kernels computed from Jupiter models 
and inversions for select regions of the interior. We use forward and inverse techniques 
developed for helioseismology. Using the example modes, we invert for sound speed and 
density perturbations. 

The three models exhibit discontinuities at interior boundary layers, and we conclude 
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that we could reasonably detect sound-speed jumps of a few percent in the inner envelope 
if the modes studied here were observed and available. For deep transitions in the core 
boundary where fewer modes probe, it is still very likely that the locations of the large dis- 
continuities could be accurately found, but the magnitudes of the jumps will be somewhat 
more difficult to ascertain because of noise. All modes have different sensitivities at different 
depths, and these results depend on a number of factors, including the parameters used in 
any inversion. While a set of simplifying assumptions was necessary to perform this analy- 
sis, these examples do provide a powerful demonstration of what to expect once confirmed 
measurements are available. 

New dedicated observations of Jovian oscillations are strongly encouraged based on these 
results, which show that most of the interesting features of Jupiter's interior are captured 
by inverse techniques. We strongly encourage observers to design a strategy to detect, at 
bare minimum, the mode set studied here. For a detailed study of the outer layers of the 
planet near the surface, even higher-degree modes that are trapped in that region would be 
necessary to measure and identify. 

We thank Patrick Gaulme for very useful comments about the paper, and acknowledge 
the suggestions by two anonymous referees that made the text significantly clearer. This 
work was supported by the NASA Outer Planets Research Program. 
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Fig. 5. — Small frequency separation 5u for each of the three models. Dashed lines are 
theoretical values (from Eq. |3]) and solid lines are computed from the model frequency 
spectrum. The angular degree pairs for consecutive orders are (1,3) and (0,2). Note the 
weak frequency dependence of the theoretical values from the term in the denominator in 
Eq.[3j 
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Fig. 6. — The lower turning points of a set of modes as a function of angular degree for each 
of the three models and for three different representative frequencies computed from Eq. HJ 
The horizontal lines deep inside Jupiter denote the boundary of the core/metallic region 
for each model, while the horizontal lines in the outer ~ 40% represent the approximate 
inner/outer envelope boundary. Associated pressures are also labeled for the outer envelope 
boundaries (model MH08 does not have one). The vertical line marks the boundary of the 
modes considered in this study (£ < 25). 
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Fig. 7. — Sensitivity kernels for squared sound speed (a) and density (b) for the three Jovian 
models. Each panel represents a different angular degree, and various values for each of the 
radial orders were chosen for illustration purposes. Also labeled is the mean eigenfrequency 
V (averaged over each model) of the corresponding eigenfunctions used to construct the three 
kernels in each panel. Note the kernels are scaled with Jupiter's radius and are plotted in 
arbitrary units. The kernels in column (a) are positive and the kernels in (b) oscillate about 
zero. 
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Fig. 8. — Example inversions f^modes with i < 25 using the SCV rnMyiaiclegfeidy the 

sound-speed jump in the outer envelope transition of that model at ~ 0.8Rj. Panels (a) 
and (b) show an inversion for a wide target function (0.2i2j), and panels (c) and (d) for 
a narrower one (0.05-Rj). On the x axis, the triangle shows the target depth ro, and the 
diamonds r ±FWHM/2 of the prescribed target function. The upper right portion of panels 
(a) and (c) provide the values of the inversion relative error (in percent of squared sound- 
speed difference) and target FWHM (in fractional radius), respectively. Panels (b) and (d) 
show the properties of the coefficients a n e(r ). The bottom plot in (b) and (d) denotes 
the relative strength of the coefficients by the size of the circles and the sign (filled circles: 
positive, open circles: negative) in £ — v space. The top plot in these panels shows the 
amplitude of the coefficients as a function of their lower turning point. 
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Fig. 10. — Example inversion u^lrjs; the LMR model with similar parameAepkaaldpeFig. [9b, 

but using a limited mode set for i < 10 and 1 < v < 3 mHz. Note the much worse agreement. 

This figure follows the style of Fig. [HJ 



-30- 



(>) 



(b) 



0.26 

0.20 



Target function 
Avg. kern, for p 

Avg. kern, for (? s 




tSJ 



^•..•'..•'%)' ,J . „ro- 



ft: 




.o* 

- • o» 



^•S- V°o«» o :oo»°°;«o»«° 





K- ' °'i~£'Z'0°'Z. -o 

•°:.-;!o-*r^>- , -.»° 



0.2 0.4 0.6 0.8 1 10 20 

Fig. 11. — Example inversions for ole^ity differences using model MH0§ to target the core 

boundary layers at about r = 0.1-Rj. The panels in column (a) show comparisons of the 

averaging kernels with the target function for varying FWHM, decreasing from top to bottom. 

The panels in column (b) show the corresponding contributions of the inversion coefficients 

in £ — v space. 
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Fig. 12. — Example inversions for dej^ity differences using model MHQ8 to target the core 
boundary layers as in Fig. [Til but varying the trade-off parameter. The panels in column 
(a) shows inversions with fixed FHWM= 0.1, but for \x decreasing from top to bottom by 
two orders of magnitude in each panel. Column (b) plots the coefficients for these inversions 
in t — v space. 



